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The structure of ''''S has been studied using delayed 7 and electron spectroscopy at GANIL. The 
decay rates of the Oj isomeric state to the 2f and Of states have been measured for the first time, 
leading to a reduced transition probability B(E2 : 2^^-0j)= 8.4(26) e fm and a monopole strength 
p^(EO : Of^Ot) = 8.7(7)xl0~^. Comparisons to shell model calculations point towards prolate- 
spherical shape coexistence and a phenomenological two level mixing model is used to extract a 
weak mixing between the two configurations. 
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'Magic' nuclei exhibit large gaps between the occu- 
pied and valence orbits. They are cornerstones of the 
nuclear structure, being used (i) to test our understand- 
ing of the nuclear forces which form these gaps and (ii) 
to model more complicated systems having many valence 
nucleons. While nuclei having 8 and 20 protons (or neu- 
trons) can be reproduced by modeling the atomic nu- 
cleus with an harmonic oscillator potential, a spin-orbit 
interaction must be added to describe heavier magic nu- 
clei. This spin-orbit interaction strongly binds nucleons 
having their angular momenta i aligned with their in- 
trinsic spin value s, denoted as £f. This leads through- 
out the chart of nuclei to regular sequence of orbits if, 
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2)i, ii, with the so-called large spin-orbit 



gaps 14, 28, 50, 82 and 126 between the lowered if or- 
bit {£=2, 3, 4, 5 and 6) and the others. Generally, in 
particular at the stability, these gaps are large enough 
to prevent excitations between occupied and valence or- 
bits and these magic nuclei are spherical. However, as 
the orbits forming the gap are separated by two units 
of angular momentum, quadrupole excitations are likely 
to develop if for some reason the shell gap is reduced. 
In this hypothesis, the development of quadrupole exci- 
tations jeopardizes the rigidity of the spherical gap and 
conduct the nucleus to deform. Consequently the doubly 
magic nuclei which have proton and neutron spin-orbit 
shell closures could become vulnerable to quadrupole ex- 
citations, as both protons and neutrons could act co- 
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herently to deform the nucleus. So far the prototypical 
deformed nucleus composed of such a double spin-orbit 
shell-closure is f|Si28 [l|. At N=28 a gradual develop- 
ment of deformation occurs between the spherical doubly 
magic 2oCa28 and the deformed f|Si28. In between these 
two extremes, i.e. in £5828, competition between spheri- 
cal and deformed shapes is expected to be present leading 
to shape coexistence 0t3| • Depending on the strength of 
the quadrupole correlations induced by the cross shell ex- 
citations either the spherical normal configuration, or the 
deformed one, becomes the ground state while the other 
configuration forms a low lying Oj state. Therefore the 
discovery and characterization of this Oj state in '*'*S28 
represent crucial information for understanding the evo- 
lution of N=28 shell gap. The non spherical nature of 
the ^''S ground state was suggested by its short /3 half- 
life and weak neutron-delayed emission probability [7| , by 
the low energy of the 2f state (1297(18) keV), and the 
enhanced reduced transition probability B(E2 : 2^— 7>0f ) 
of 63(18) e^fm* [8]. However the 2+ and B(E2) val- 
ues are intermediate between a rigid rotor and a spher- 
ical nucleus. It suggests a possible mixing of spherical 
and deformed shapes which can be deduced by studying 
the properties of the Oj isomer at 1365(1) keV observed 
in [9]. Already the study of a 7/2" isomer in ^^S p, Q 
has shed light on shape coexistence in the N«28 region. 
Other cases of shape coexistence around shell closures 
have been reported in [l^, Ull ■ 

The present letter reports on the determination of the 
monopole strength p^(EO : Oj— >-0f) and the reduced 
transition probability B(E2 : 2^^— >-0j) in *^S, extracted 
from the measurement of the half-life and the branch- 
ing ratio between the EO and E2 decay of the isomeric 



Oj state. These pieces of information were obtained by 
using combined 7 and electron delayed-spectroscopy and 
are used to demonstrate the shape coexistence in '^^S. 

The experiment was carried out at the Grand 
Accelerateur National d'lons Lourds (ganil) facility. A 
primary beam of "^^Ca at 60A-MeV (I^2e/iA) impinged 
onto a 138 mg/cm^ Be target to produce neutron-rich 
fragments. They were separated by the lise3 spectrom- 
eter [iJI using an achromatic 100 mg/cm^ Be degrader. 
The magnetic rigidity was set to optimize the transmis- 
sion of the "'^S nuclei, produced at a rate of 200 sec^^, 
with a momentum acceptance of ±1.45%. Fragments 
were identified on an event by event basis by means of 
their energy loss and magnetic rigidity (Bp) values. The 
Bp was obtained from the position of the fragments at 
the dispersive focal plane given by a multi-wire propor- 
tional chamber (caviar) [IJ]. The selected nuclei were 
implanted in a 125 /im kapton foil tilted at 20 degrees 
with respect to the beam axis. Before the foil, a stack of 
Si detectors, including a position-sensitive one, was used 
to adjust the implantation depth and to reconstruct the 
position of the ions in a plane perpendicular to the beam 
axis. A thick Si detector located downstream to the im- 
plantation foil was used as veto. The 7 and electron decay 
events were registered up to 20 /isec after the implanta- 
tion. Electrons were detected in four cooled 45*45 mm^, 
4 mm thick Si(Li) detectors, placed 20 mm above and 
below the beam axis. The 7-rays were measured by two 
clover Ge detectors of the exogam array located on the 
side of the implantation foil, at a distance of 25 mm to 
the beam axis. The use of a parallel beam optics along 
2 meters length enables to derive the ion implantation 
profile on the kapton foil from the position-sensitive Si 
detector. This ion profile, the geometry of the detec- 
tors and that of the chamber were used as ingredients 
in a GEANT4 simulation to derive the electron (eg- ) and 
7 (e-y) efficiencies. The simulated efficiencies compare 
well with the ones obtained with calibrated sources of 
^°^Bi and ^^^Eu placed in calibration runs at 6 difi^erent 
positions on the implantation foil. Using these compar- 
isons, 6^=3.06(5)% and eg- =13.3(6)% were adopted for 
a gamma-ray of 1329 keV and an electron of 1362.5 keV, 
respectively 
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The decay of the J to the ground-state (EO) proceeds 
through the emission of an internal conversion electron 
(IC) and by internal pair formation (IFF). The electron 
spectrum, following the implantation of a ^"^S nucleus, is 
shown in Fig. [T] A single peak is observed at 1362.5(10) 
keV corresponding to an excitation energy of 1365(1) keV 
for the O2 state, after having corrected for the binding en- 
ergy of the K electrons in the ''^'^S nucleus. The integral of 
the peak is I^- {E0)= 148(8) xlO^. The low energy part 
of the spectrum is well accounted for by the pair forma- 
tion (IPF) mechanism in which electrons and positrons 
share an energy of 1365-(2*511)=343 keV. The fit of the 
electron time distribution (insert of Fig. [1]) leads to an 
half-life of 2.619(26) /isec, which agrees with the value of 
2.3(5) /isec reported in [9i]. 
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FIG. 1: Electron energy spectrum obtained from the Si(Li) 
detectors following the implantation of ''''S nuclei. The peak 
at 1362.5(1.0) keV corresponds to the 0;j"— ^0^ EO transition. 
The low energy part is due to pair creation. Insert: Time 
distribution of the 1362.5 keV electron transition from which 
a half-life of 2.619(26) /isec is extracted. 



The 0^^2]^ decay branch, occurs through a strongly 
converted E2 transition at 36(1) keV, an energy below 
the experimental threshold of the detection system. The 
energy of this unobserved transition is derived from the 
measured energies of the J and 2^ states, the latter be- 
ing obtained from the observation of a delayed 2^— i>0^ 
transition at 1329.0(5) keV (half-life of 2.66(23) psec, 
in agreement with the value obtained from the elec- 
tron spectrum) which follows the Oj— ?>2j'" decay. The 
1329 keV energy agrees with the value of 1297(18) found 
in Ref. fs"]. The yield of the Oj-J'2:;'" transition, /^(£'2), 
has been extracted from the number of delayed 2^— J'O^ 
7-rays. As can be seen in the insert of Fig. [51 this tran- 
sition is contaminated by the 1332.5 keV 7-ray of ^"Co 
arising from the activation of the last selection slits of the 
spectrometer, which also produce a 1173 keV 7-ray with 
the same intensity. The number of counts in the 1329 
keV peak has been obtained by fitting the 7 spectrum 
with two gaussians, the intensity of the 1332.5 keV tran- 
sition being constrained by that of the 1173 keV 7-ray. 
The resulting Ij{E2) is 56(3) xlO^ 

The decay of the Oj state occurs through E2 and EO 
transitions, the ratio of which is expressed as : 
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In this expression, the electronic factors for pair for- 
mation and internal conversion have been extrapolated 
for a nucleus with A=44 from Ref. [15l - ll7l | to be 
r2^P^=1.495xl0^sec-i and f}/c=1.1125xl0^sec-\ re- 
spectively. A value of 3.6xl0~^ has been taken for the 
conversion coefficient aco™(2^^0^) [l^|. Using the ex- 
perimental values of electron I^- (EO) and 7-ray I-y {E2) 
yields derived above, the resulting branching ratio is 
R=0. 163(13). The p^{EO : Oj^Of ) and B(E2 : 2+^Oj) 
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FIG. 2: Part of the delayed gamma energy spectrum following 
the implantation of '''*S nuclei. Peaks from the /3 decay of ^^K 
(1158 keV) and ^°Co (1173 and 1332.5 keV) are identified, the 
latter overlapping with the 1329 keV 2^— s-O^ transition of ''^S. 
The deconvolution of this doublet is shown in the insert. 



values are obtained using the follovifing equations 
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f|S28 using the antoine code [19[ and the up-to-date 
SDPF-u interaction that accounts remarkably well for nu- 
clear structure in this mass region [20|. The full sd (fp) 
valence space has been considered for protons (neutrons) 
using standard effective charges e^ = 1.35 e (e,y — 0.35 e). 
The results gathered in Table H] show a good agree- 
ment with the experimental values, the only exception is 
a somewhat larger calculated B(E2 : 2]^^0j) value than 
measured. Nevertheless, both experiment and calcula- 
tion agree with the fact that the 2'^ state connects much 
strongly with the 0^ state than with the 0^ one. In- 
deed, the experimental B(E2 : 2+^05^)/B(E2 : 2+^0^) 
ratio is 7.5 whereas the calculated one is 3.2. Calculated 
excited states connected to these two 0"*" states are pre- 
sented in Fig. 13] with their intrinsic quadrupole moments 



TABLE I: Experimental and shell model values for the exci- 
tation energies, in MeV, and reduced transition probabilities 
B(E2), ine^fm*, of ''''S. 
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exp. 1.329(1) 1.365(1) 2.335(39) 

SM 1.172 1.137 2.140 
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Using the measured branching ratio R, the half- 
life value Ti/2(0j) and aeo™(0j->2+)==10.94(l) ex- 
trapolated from Ref. [1^, the monopole strength 
p^(E0 : 0;^— >0]^) and the reduced transition probability 
B(E2 : 2+"'-^0j) have been determined to be 8.7(7) xlO^^ 
and 8.4(26) e^fm^, respectively. 

The values of E(Oj) =1365(1) keV and 
p^(E0)=8.7(7)xl0-'^ are the smallest measured in 
this mass region, pointing to a weak mixing between 
the 0^ ground state and the 0^ isomer and therefore Q 
to shape coexistence. In case of a large mixing, these 
states would repel each other to exhibit a large energy 
spacing and a larger p^(EO) value. To obtain further 
understanding on the nature of the shape coexistence, 
data are compared to shell model calculations. 

Shell model (SM) calculations have been performed for 



FIG. 3: **S level scheme calculated within the present SM 
approach (left), compared with available experimental data 
(right). E2 transition probabilities (in e^fm*) are reported 
on top of black arrows and intrinsic quadrupole moments (in 
efm?) are shown in light gray on the right side of calculated 
levels. The ground state of the nucleus is head of a rota- 
tional band (/3 ~ 0.25) and coexists with the rather spherical 
low-lying 0^ isomer. Calculated values of the N=28 gap and 
correlation energies (in MeV) are given for even-even N=28 
isotones. 



Q. For sake of clarity only the states of present interest 

are shown in this picture. Remarkable is the presence 
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states on top of the 0^ ground state 
connected by large B(E2) values. These states present 
equal Qo values of about 60 efm?. These two features 
characterize the presence of a rotational band from an 
axially deformed nucleus with /3 ~ 0.25. The 2^ state 
at 2.14 MeV has a smaller intrinsic quadrupole moment 
Qq—-0.3 efm? compatible with a spherical shape. A can- 
didate for the 2^ state is proposed at 2335(39) keV by 
placing the previously reported 988 keV transition [2l| 
on top of the 0^ or 2'1 state. Hence SM calculations 
suggest a prolate-spherical shape coexistence in ''^S. 

A detailed analysis of the components contributing to 
the total energy of the 0+ states has been performed in 
order to deepen our understanding on the evolution of the 
collectivity from 2oCa to 'WSi. Within the SM framework, 
the total Hamiltonian can be separated into its monopole 
(i.e. spherical mean-field contribution) and multipole 
(i.e. correlations mainly of pairing and quadrupole type) 
parts |22l |. As can be seen from the values reported 
in Fig. 131 correlations strongly increase from the dou- 



bly magic '*®Ca (~2 MeV) down to the exotic deformed 
''^Si (~18 MeV), while the size of the N=28 sheh gap 
gets shghtly reduced [23] ■ This increase of correlations 
is favored on one hand by neutron quadrupole excita- 
tions across the N=28 gap between the /7/2 and P3/2 
orbits 123], and on the other hand, by the degeneracy of 
proton Si /2 and ^3/2 orbits and excitations from the d^i2 
shell [ij, l24-26J. In both cases, quadrupole correlations 
are favored by the fact that occupied and valence states 
are separated by two units of angular momentum. With- 
out considering multipole contributions to the OJ" and Oj 
states in '*^S, both levels are found to be quasi-degenerate 
in energy, and the ground state of ^"^S is spherical. A 
gain of 1.5 MeV from the multipole energy brings the de- 
formed configuration at the minimum of binding energy, 
while the spherical configuration corresponds to the ex- 
cited state. Similar multipole effects energetically favor 
the oblate O"*" state in ^^Si which is predicted to coexist 
with a prolate 0+ state [13] at 1293 keV. 

The shell model calculation uses an Harmonic Oscil- 
lator basis for the description of the atomic nucleus. 
From the definition of the EO operator, the calculated 
EO transition between states of the same harmonic oscil- 
lator shells, as for protons in the sd shells and neutron in 
the fp shells, is strictly zero. Therefore, in order to shed 
light on the amount of mixing between the OJ'2 states 
and to deduce their shape before mixing, we use a phe- 
nomcnological two interacting levels model. We assume 
two spherical (S) and deformed (D) states before mixing 
which interact to produce 0^ and Oj states defined as : 



]0i) =cos6I]0d) +sine\Qs) 
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(4) 
(5) 



where 9 is the mixing angle. The E2 transition be- 
tween the 2^ and Oj (or Q^) states being mainly due 
to the D-component of these 0''' states, it follows that 
B(E2 : 2+^0j)/B(E2 : 2+^Q+)r^ta.iY^{e) (eq. 2 of [13]). 
A mixing amplitude tan'^(6')=0.13 is deduced from the 



experimental B(E2) values whereas the shell model gives 
a somewhat larger value of 0.24, both being smaller than 
the case of a maximum mixing (tan^(0)=l). Therefore, 
the shape coexistence is found to be more pronounced 
experimentally than calculated by the SM. The magni- 
tude of the monopole matrix element can be written as 
a function of the mixing amplitude and of the difference 
of shapes, Ps and /?£i, between the two configurations 
before mixing ^, p'^{EQ)={2,Ze/'i'iTfsin'^6cos'^e{l3jj - 
/3g)^. Using the experimental mixing amplitude value 
of iSiT?{9) in this equation, the experimental monopole 
strength is reproduced only when deformations /3zj ~ 
0.274 and /3s=0 are assumed. The deformation param- 
eter Pd is in close agreement with the values obtained 
after mixing from Coulomb excitation experiment [8|, 
/3 = 0.258(36), and from the shell model calculations, 
/3=0.25. Altogether these values again point towards a 
deformed-spherical shape coexistence in ^^S. 

To summarize, electron and 7 delayed-spectroscopies 
have been used to determine the monopole strength 
p2(E0 : 0j^-0;i^)=8.7(7)xl0"3 and the reduced transi- 
tion probability B(E2 : 2+^0j)=8.4(26) e^fm^ in the 
16^28 nucleus. Using these values, the earlier measured 
13 (E2 : 2'^^Q'^) Q, shell model calculations and a two 
level mixing model, it is found that ^^S exhibit a shape 
coexistence between a prolate ground state (/3 ~ 0.25) 
and a rather spherical Oj state. This establishes how 
the onset of collectivity progressively develops between 
the spherical ^Ga. and the deformed f|Si nuclei. This 
study completes uniquely the understanding of the shell- 
breaking mechanism at the spin-orbit closed-shell N=28, 
which is as well of importance for the evolution of other 
shell gaps having the same origin. 
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